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Ph' Abstract 

The evolution equations for Green functions of the Reggeon Field Theory (RET) are 
I equivalent to those of the inclusive distributions for the reaction-diffusion system of classical 

particles. We use this equivalence to obtain numerically Green functions and amplitudes of 
CN , the RET with all loop contributions included. The numerical realization of the approach is 

described and some important applications including total and elastic proton-proton cross 
' sections are studied. It is shown that the loop diagram contribution is essential but can be 

. imitated in the eikonal cross section description by changing the Pomeron intercept. A role 

of the quartic Pomeron coupling which is an inherent part of the stochastic model is shown 
to be negligible for available energies. 
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1 Introduction. 

The Reggeon Field Theory (RFT) [1] is successfully used up to nowadays for describing 
the soft part of the strong interaction dynamics at high energies. It is based on the 
very general principles which should obey s-channel amplitudes such as analiticity, 
crossing and t-channel unitarity. Particular field theories like QCD should obey the 
same principles. The Regge-like behavior was observed in QCD and some of the 
RFT components which are not deducible from the first principles like the Pomeron 
intercept were obtained within certain approximations. These approximations of the 
QCD however are questionable to apply to the processes which do not involve high 
momentum transfer. Moreover, it is rather complicated to fulfil the reqirements of full 
many-particle t-channel unitarity which is very important for high-energy dynamics. 
Due to the generality and phenomenological success of the RFT in many aspects it can 
be viewed as complementary to the QCD and should be treated on the equal grounds 
with the QCD-induced models. 
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For the description of the experimental data in the reggeon approach one usually 
uses the simplest reggeon diagrams corresponding to the pomeron pole and a set of 
nonenhanced regge-cut graphs providing s-channel unitarity. However for the case of 
the supercritical pomeron (with the intercept larger than one) the interaction between 
pomerons may be essential and can become especially important at the LHC energies 
so the contribution of loop reggeon diagrams should be estimated. In this paper we 
present such estimates for the RFT with triple and quartic pomeron vertices using the 
mathematical equivalence of this theory to dynamics of stochastic system of classical 
particles [2]. 

The Lagrangian density of the Reggeon Field Theory is formulated by using the 
non-relativistic Pomeron field </>(?/, b) 

^ = \<P\% - t)<P - «'(Vb0^)(Vb0) + A0V + Ani , (1) 

where the rapidity y plays a role of time, b represents position in the impact parameter 
plane, parameters a' and A describe the linear Pomeron trajectory, ap{t) = (1 + A) + 
a't, and L^t contains the interaction terms. The usual choice for the interaction terms 
is the option when it generates only triple Pomeron vertices, however terms of higher 
order in and (p^ can also enter. 

It was shown in [2] and discussed in the subsequent papers [3, 4] that the classical 
stochastic system of particles (called "partons" for brevity^) which undergo diffusion 
in 2-dimensional space and interact in a specific way (the particles can split, die, fuse 
or annihilate in pairs) allows a field-theoretical realization corresponding to the RFT 
Lagrangian (1) with 

Cint = I r3p0V(0^ + 0) + + Jo0^ + Jy0 • (2) 

The interaction includes both triple (nonhermitian) and quartic terms, the source 
term jo(0,b) gives the initial condition for the ^/-evolution ("parton" distribution in 
the projectile) and the term jY{Y, b) describes the partonic structure of the target and 
is necessary for calculating the scattering amplitude at total rapidity Y . The values of 
the Reggeon couplings in (2) are connected with the rate parameters of the stochastic 
model. Equality of couplings r.^p for both triple terms in (2) imposed by Lorentz 
invariance requires tuning of the parameters. The multiparticle inclusive distributions 
of the stochastic system turn out to be identical to the Green functions of the RFT. 
This correspondence may be used to compute the amplitudes of that specific model of 
the Reggeon field theory with the account of all loops. 

In this paper we present a numerical realization of the stochastic approach and 
give some estimates of the loop diagram contribution for the RFT in ti = and d = 2 
dimensions. In the next section we recall the notations and the definitions of the 
stochastic approach and outline the correspondence between the parameters of the 
RFT and that of the stochastic approach. Section 3 is devoted to the description of 
the numerical method used for the computations of the amplitude. In Section 4 we 
analyze asymptotic behaviour of the Pomeron Green function in d = and d = 2 
dimensions for various sets of parameters. The role of the quartic Pomeron coupling 

"'^We note that these classical partons should not be confused with the partons in the Fock wave function or the 
color dipoles in the QCD induced cascade models [5, 6]. 
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is discussed. In Section 5 the correspondence between the parameters of the stochas- 
tic model and the couphngs of the Reggeon field theory is established by comparing 
contributions of simplest diagrams. We also discuss the restrictions which one must 
impose on the numerical algorithm to reproduce the contribution of the limited sets 
of the Reggeon diagrams which enter in the (quasi)eikonal approximation. Then the 
calculated contribution to proton-proton cross sections of the full set of the Reggeon 
diagrams is compared to the standard calculations in the quasieikonal approximation 
and the role of the Pomeron interactions is estimated. The effective change of the 
Pomeron intercept due to account of these interactions is found. Concluding remarks 
are presented in the last section. 



2 The stochastic model 

Here we briefly describe the stochastic model of the RFT (also known as reaction- 
diffusion approach) as it is formulated in [3] . We follow the evolution of the system of 
classic particles, "partons", which interact in a certain way. The partons are allowed 
to move randomly in the transverse plane, and this diffusive motion is described by a 
diffusion coefficient D. The particles can split, A ^ A + with a splitting probability 
per unit time A, or die, A — !■ 0, with a death probability mi. When two partons 
are brought within the reaction range due to the diffusion, they can pairwise fuse, 
A + A ^ A, oi annihilate, A + A ^ ^. To describe the rate of these processes we use 
dimensional fusion and annihilation constants, ay and ami'- 

^^ = j d%py{h) , am2 = j d^bpm^ib) ■ (3) 

The dimension of and a^, thus is length squared. The characteristic scale, that is 
the effective interaction distance, is determined by the functions Pu{b) and pmiip). It is 
assumed to be small compared to the characteristic scale of the problem (e.g. proton 
radius).^ 

Throughout the evolution the system of partons is described by a set of symmetrized 
probability densities pjv(l/; bi, . . . , b^r) of finding exactly N partons at evolution time 
y at positions {bi, . . . , b^v} in the transverse plane. For brevity we use the notation 
Bm for the set {bi, . . . , bjv} and dB]^ for dhi . . . dhjq. The probability densities are 
normalized according to 

■^jdBt,pt,{y;Bt,)=pt,{y), (4) 
where PN{y) is the probability of having exactly partons at the evolution time y, 

The important quantities which have a direct analog in the Reggeon field theory 
are the inclusive s-parton distributions, which correspond to observing s partons at 
positions Zg = {zi, . . . , z^} and integrating out all the rest. They are defined as 

fs{y;Zs) = X] _ / dBN PN{y;BN)YlS{zi-hi). (5) 

N ^ ^''-^ 1=1 



^ Note that in ref . [3] the fusion and annihilation constants (denoted there as v and m2 respectively) are defined 
with respect to a single parton and thus are twice less compared to Ci/ and (Jm2- 
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Due to the symmetry properties of pNiy] ^n) the inclusive distributions fsiy] Zg) are 
also symmetric functions of their variables. As follows from the normalization of 
the the inclusive distributions are normalized to the factorial moments /ig of the 
distribution PN{y)- 

f ^ N\ 

As it is pointed out in [3] the evolution equations in time y for the inclusive distri- 
butions fs{y; Zs) 

^fs{y; Zs) = D Vlfsiy; ^.) + (A - m,)sfs{y; Z,)+ 

+xJ2 fs-i{y, zi'^)6{zk - zi)-uJ2 fs{y; zM^k - z,)- 

k,l=l k,l=l 

s 

-(2m2 + z/)J]/.+i(y;2f),Zfc,Zfc) , (s = l,2,...) (7) 

k=l 

exactly reproduce the evolution equations in rapidity for the Green functions of the 
RFT with the interaction term (2). This means that the inclusive distributions them- 
selves are proportional to the Green functions of the RFT in convolution with the 
particle-m-Pomeron vertices of the RFT N'mi'^m)'- 

fs{y]Zs) (X^ / dX,nJ^m{'^m)Gmn{0]Xm\y]Zn) ■ (8) 

m 

where the Green functions Gmn{0] '^m\y] Zn) describe the transition of m Pomerons 
at rapidity y = into n Pomerons at rapidity y. The proportionality factor takes 
into account the dimension of the inclusive distributions and should be specified at 
calculating Regge amplitudes (see below). The vertices ^fm,{,'^m) determine the initial 
distribution of partons in number and transverse positions at zero time y = which 
are used as a starting point for the evolution of fs- We will attend to this point in 
section 5. 

The amplitude of the RFT as a function of the total rapidity Y and impact param- 
eter b is expressed in terms of the set of the inclusive partonic distributions fs{y, Zg) 
for the projectile A and fs{y — y, Zs) for the target A taken at arbitrary intermediate 
rapidity y through the procedure of linkage: 

T{Y) = {A\T\A) = ^ / dZsdZsfsiy, Zs)fs{Y - y- Z^) \{ g{z, - z, - b). 

s=l i=l 

(9) 

This procedure realizes a linkage of partons (Pomerons) of the target and the projectile 
through a convolution with the linking functions g{. . .). These functions are narrow 
functions of the 5- function type^ normalized to some parton size e = vra^. They will 
be specified in section 3. 



The original paper [3] had J|(5'^^(zi — Zi — b) in (9) instead of more general expression n5(-^s' ~ Zi — h), 
however this restriction is not essential and has mainly the computational significance. The only important thing 
is that the characteristic scale e of the functions g must be small compared to the scale of the problem. It differs, 
in general, from scales of the functions pi, and ■ 
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The formula (9) links two partonic cascades coming from the target and the pro- 
jectile at the rapidity value y. It was proved in [3] that, provided the scales of the 
functions g{h), Pu{b) and Pm2{b) are small, the convolution (9) depends only on the 
overall rapidity/evolution time Y but not on the position of the linking point y ii a 
certain condition on the model parameters is satisfied. Namely, 

A J g{h)dh = j p,nMdb pMb , (10) 

or, 

1 

Ae = 0-^2 + -CTu ■ 

This condition in fact corresponds to the equality of the Pomeron fusion and splitting 
vertex. 

The relation between the Reggeon field theory parameters and the ones of the 
stochastic model is given in table 1. 



Table 1: Relation between the parameters of the RFT and those of the stochastic approach. 



RFT 


stochastic model 


Rapidity y 
Slope a' 
A = a{0) - 1 
Splitting vertex r^p 
Fusion vertex r^p 
Quart ic coupling x 


Evolution time y 
Diffusion coefficient D 
X — mi 
Av^ 

(0-m2 + \<^u)/\^ 



Since the inclusive distributions (5) are in direct correspondence with the exact 
Green functions of the RFT with the account of all loops (the model characteristic scale 
of py and Pm-2. plays a role of the cutoff for the RFT Green functions), the equation (9) 
can be used for obtaining the exact numerical answer for the RFT scattering amplitude. 



3 Numerical method. 

The inclusive distributions js and js which enter (9) are defined by the evolution of the 
partonic subsystems with a certain initial distribution of partons in their number and 
transverse positions. Basically one may assume a straightforward way of calculating 
the interaction amplitude according to (9) by making a Monte-Carlo average over the 
initial conditions and computing /s and explicitly. This method, however, seems to 
be numerically unaffordable. For example, besides the evolution of partonic system and 
computing the distributions on some grid which one has to introduce, it is necessary 
to perform symmetrization of the /^'s which (provided the grid has fairly small cells) 
is extremely time consuming, or, alternatively, increase the number of Monte-Carlo 
simulations in the same proportion. The corresponding number of operations can be 
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estimated as M!/(M — s)!, which is the number of choices of s out of M points of the 
grid together with number of transpositions. 

Let's rewrite (9) exphcitly indicating the averaging over the initial distribution of 
partons: 

s=l n,n 

I dZJZJs{y-,Bn\Zs)fsiY-y-Bf,\Zs)l[g{z,-i,-h) (11) 

Here fsijj] Bn\Zs) and fs(Y—y; Bn\Zs) stand for the s-parton inclusive distributions ob- 
tained upon the evolution of the initial sets of projectile-associated partons at positions 
Bn and of target-associated partons at Bn. The straightforward way of calculations 
described above consists in averaging over the sets Bn and Bn prior to the integration 
over Zs and Zg. 

Another way of calculations can be realized making averaging over {ra, Bn} and 
{n, Bn\ after the convolution in Zg and Zg. In other words, the Monte-Carlo aver- 
aging is done for the amplitude computed on the "samples", particular realizations 
of the probability distributions which correspond to the set of inclusive distributions 
fs and fg. Each sample is generated by evolving a set of partons in rapidity starting 
from the random number and initial positions of partons generated according to the 
distributions P{n,Bn) and P{h,Bn). 

A set of N partons with fixed coordinates = {xi, . . . ,XAr} in the transverse 
plane corresponds to the following set of the s-parton inclusive distribution functions 

fgiZg) is = i...Ny. 

fg{Zg)= Yl 5(zi-x,J...5(z,-x,J. (12) 

{x,^,..,x,^}eA'jv 

The sum goes over all the ordered subsets of s variables out of which ensures the 
normalization of fg and its symmetry properties. The same is for fg{Zg) where the 
parton coordinates jji G 3^^. After substituting (12) the sample amplitude (cf. with 
(11)) reads: 

N 

Yl 9{K - yn)---9{^is - Yjs)- (13) 

jl,...,js = l 

Here A^ and A^ are the numbers of partons in the samples which come from the evolution 
of the projectile- and the target-associated sets of partons and Nmin = min(A^, A^) is 
the minimal number of partons in these two samples. 
The values 

9ij =g{^i- Yj) (14) 

computed on the samples of partons form a matrix which can be looked upon as the 
adjacency matrix of the bipartite graph with weights gij. 



^sample = ^ ^ 

s=l ii,...,is = l 
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For a given s in the sum (13) there are s! identical terms which correspond to 
permutations of the pairs {xj,yj}. This is since the summation in (13) is goes not 
only over the subsets {ii, . . . ,is} G {1, • • • , N} and {ji, . . . ,js} G {1, . . . , A^} but also 
over the permutations within these subsets due to symmetrized form of fs{Zs). Taking 
this into account we arrive at the expression 

^ min 

^sample = ^ ("l)""^ ^^1^1 " " " ^'^^^^ ^^^^ 

s=l ii<i2...<is ji,---,js 

It means that for each term of the sum with fixed s one has to take all possible 
products of s elements of the matrix G = \\gik\\ so that the elements from the same 
line or column would not enter same product. The products with the odd number of 
multipliers enter with +1, the ones with the even - with —1. Speaking in terms of 
mathematical definitions, this is a sum of the permanents of all the square submatrices 
of the matrix G = {gij} with alternating sign. The submatrices are obtained from the 
initial matrix G by removing the relevant number of lines and columns. 

The permanent of a matrix (see e.g. [7]) has a definition analogous to that of the 
determinant, but without the alternation of the sign in the sum: 

n 

Per{A}= J2U^^^ (16) 

Despite its similarity with the determinant, the permanent is quite an inconvenient 
construction for a numerical computation. The fastest known algorithm is given by 
the Ryser formula: 

Per(A) = (-l)'^ J2 (-l)'"'nE«^^' (17) 

5C{l,...,n} 1=1 j&S 

where summation goes over all the subsets S* of {1 . . .n} and 15*1 denotes the number 
of elements in the subset. The number of operations for computing according to this 
formula is around 0(2"n). This is much less than 0(n!) needed for the straightfor- 
ward computation according to (16), however still incomparable with O(n^) operations 
needed for computing the determinant. 

One can easily generalize Ryser formula for the case of the amplitude Tsampie in 
(15): 

siC{l,...,7V} S2C{1,...,7V}, «Gsi \j£s2 J 

\si\<\s\\ 

The estimated number of operations here is around 0(n4") for a sample, which is 
much less than the number of operations required for symmetrization of the full set of 
inclusive distributions fs{2s) on the grid in the straightforward way of the amplitude 
calculation. 

Calculations based on the Monte-Carlo averaging of (18) allow to to set an arbi- 
trary rapidity for the linkage point since the amplitude T depends only on the overall 
rapidity, that is the sum of the evolution times of two samples, as proven in [3]. The 
freedom in the choice of the linkage point can be quite convenient for calculating the 
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amplitude at large rapidities when the Pomeron intercept exceeds unity. In this case 
the number of partons in the set grows exponentially with the evolution time and at 
some point the Monte-Carlo evolution of the set can become very numerically expen- 
sive. So choosing the linkage point e.g. as a half of the overall rapidity can in principle 
compensate the difficulties connected with the computing of the permanents. Besides 
this, (18) also serves as an expansion in the number of Pomeron exchanges between 
the upper and lower blocks of the amplitude at the rapidity of the linkage point. 

If the application which we are using the stochastic approach for do not involve this 
expansion and the overall rapidity is not too large, a much less numerically expensive 
method can be used. Taking advantage of the independence of the amplitude (9) on 
the linkage point we use the fact that for zero evolution time the inclusive distributions 
fs{y = 0,^s) are proportional to the particle-mPomerons vertices (cf. (8)) and hence 
are known. Substituting (12) into the (9) together with the known functions fs{y = 
0,Zs) and as before taking into account the symmetry properties of fs{Y,Zs) in (12) 
we have: 

TV 

Tsample(b) = ^(-l)^^V.e^ J2 p(x,, - b, . . . , X,^ - b) . (19) 

s = l ii<i2...<is 

Here jig denotes the factorial moments of the distribution which corresponds to the set 
fs{y = 0, Zg) and Ps stands for its coordinate dependent part normalized to unity. In 
writing (19) we also assumed that the functions (^(x) are narrow and replaced them 
by e5(x). In case of no correlations in the vertex which for example is the case for the 
eikonal vertices ps{Zs) = YlPAi'^i) (19) reads 

N 

Tsampic(b) = ^(-l)^-V.e^ ^ p^(x,, -b)...p^(x,^-b). (20) 

s=l ii<i2...<is 

The amplitude is obtained as the Monte-Carlo average of (19) or (20) over the sets of 
partons. 

For the numerical calculations one must also fix the form of the distributions (b) 
and Pu{b) together with making a specific choice for the functions g{b) in the convolu- 
tion (9). We take these distributions in the following form: 

Pm2(b) = m2 Q{am2 - |b|); (21) 
p,{h) = z/e(a,-|b|); (22) 
g{h) = A;e(a-|b|); (23) 

where we have introduced explicitly the interaction and convolution scales; the con- 
stants u and m2 play a role of fusion and annihilation probabilities for a pair of partons 
per unit time provided they are close enough in the transverse plane; k is some ap- 
propriate numerical factor ^ . The boost invariance condition (10) in these notations 
reads: 

Xkna^ = m27ra.^2 + |7ra^ (24) 
It allows to vary the scale a preserving normalization of the function g, J (Pbg{h) ~ k tto^ ~ e. 
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while triple and quartic coupling constants are 




(25) 



For a given value of the triple coupling r^p by taking a higher value for the parameter k 
one can release the scale a and set it small enough compared to the characteristic scale 
of the problem, which is, as mentioned, the condition of applicability of the stochastic 
approach. Though three radii can in principle differ, in our numerical calculations 
for clarity we take all three the To reduce the number of free 

parameters in the model we also choose k = 1. The boost invariance condition for our 
choice takes the most simple form: 



The elastic amplitude is computed as a Monte-Carlo average of the expressions 
(13) or (19) over the sets of partons which underwent stochastic evolution over the 
given time (rapidity) interval. The way of preparing these sets is straightforward. 
First, the number of partons and their position in the transverse plane are generated 
acccording to the distribution defined by the particle-m-Pomeron vertices. After this, 
based on the position of partons in the transverse plane we define a probabilities 7j 
of splitting, death, fusion and annihillation of partons per unit time provided they 
preserve their positions. The survival probability (probability of having no change in 
the number of partons) as a function of time y in this case is a simple exponential 
p{y) = exp(-^7i?/). 

We generate time interval Sy before the change in the parton number according 
to the distribution p{y). If the time interval does not exceed yo = O.la/y/D (10% of 
the characteristic diffusion timescale) we define the type of the changing number event 
(toss it according to the total splitting, death, annihillation and fusion probabilities) 
and a candidate parton or a parton pair. Then we randomly change the positions of all 
the partons in the set according to the diffusion law with the time Sy, do the necessary 
change in the number of partons and recalculate the number change probabilities 7j 
with the new parton number and the partons new coordinates. 

If the time interval 6y exceeds yo we change the positions of partons in the set 
according to the diffusion law with the time yo and recalculate the probabilities 7^. 

The steps are repeated until we reach the desired overall evolution time. 

4 The Pomeron propagator 

Now it is straightforward to attend to the numerical studies of the RFT. As a first 
example we compare exact numerical results for the Pomeron propagator in zero- 
dimensional RFT and RFT with the account of diffusion. The zero- dimensional case 
has been studied before both analytically in a series of works starting from [8] and 
numerically ( e.g. [9]) for zero value of the quartic coupling. The most important 
result was the decrease of the two-point Green function or the Pomeron propagator 
(describing transition of one Pomeron into one Pomeron) with rapidity although the 
perturbative propagator has an exponential growth. It has also been reported in [4] 



A = 7712 + 1^/2 . 



(26) 
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that in the zero- dimensional theory with the 2—7-2 couphng a constant asymptotic 
behaviour of the propagator is possible for some special choice of the quartic coupling, 
namely x = '"sp/^ our notations. 

In terms of the stochastic model the propagator is the single particle inclusive 
distribution /i(?/,b) for the state which have evolved out of a single parton at y = 
(see (8)). The transition to zero-dimensional case is made in our numerical approach 
by setting the diffusion coefficient equal to zero. Within our calculation scheme we 
investigate the role of the quartic Pomeron coupling in the asymptotic behaviour of 
the propagator or amplitude in both zero and two transverse dimensions Reggeon field 
theory. 

We start with the zero dimensions case. Our aim is to illustrate the role of the 
enhanced diagrams simply comparing the simple pole contribution to the propagator 
with the result of the full numerical computation. 

Here the propagator as a function of rapidity is simply the average number of 
partons as a function of the evolution time. The evolution starts with a single parton. 
In fig.l we present calculations for different sets of the parameters A, z/, mi, m2 (see 
table2). The values of the parameters correspond to the same value of the intercept 
A = A — mi = 0.1, satisfy the Lorenz invariance condition in zero dimensions A = m2 + 
Z//2 (cf. eq.(66) of [3]) and differ by values of the couplings r^p = A and x = (m2 + i^)/2 
only. 

Table 2: Parameters of the model and the RFT couplings 



Set 


A 


V 


mi 


m2 


A 




X 


1 


0.1 


0.2 








0.1 


0.1 


0.1 


2 


0.1 


0.1 





0.05 


0.1 


0.1 


0.075 


3 


0.1 








0.1 


0.1 


0.1 


0.05 


4 


0.15 


0.3 


0.05 





0.1 


0.15 


0.15 


5 


0.15 





0.05 


0.15 


0.1 


0.15 


0.075 



As one can see from fig. 1, the asymptotic behaviour of the propagator for zero- 
dimensional theory strongly depends on the relative values of m2 and v or, in terms 
of the RFT on the relation between the intercept and the couplings. The only possi- 
bility to have constant asymptotic behaviour distinct from zero within the stochastic 
approach is to put m2 = mi = as in the set 1. This corresponds to equal triple 
and quartic coupling with A = A (where A -|- 1 is the Pomeron intercept) within zero 
dimensional RFT which evidently satisfies the condition necessary for the constant 
propagator asymptotic behaviour outlined in [4]. The constant asymptotical value 
corresponds to the eq. (27) of [3] /ii(oo) = e/(e — 1) ~ 1.582. The same result is 
obtained by explicitly solving numerically the differential equation (eq. (6) of [3]) for 
generating function of the probability distribution. 

Let us emphasize that all parameters of the stochastic model should be positive. 
Another strong restriction is given by the Lorentz invariance condition which excludes 
some increasing regimes of the stochastic model. 

For the two dimensional case we take the same values of the stochastic model pa- 
rameters as for the zero-dimensional one. Additionally we introduce the partonic inter- 
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Figure 1: The full Pomeron propagator in zero dimensions as a function of rapidity. The stochastic 
model parameters corresponding to the sets 1-5 are in the text. 

action distance a = 0.05 fm and the diffusion coefficient a' = 0.01 fm^ = 0.26 GeV~^. 
The Pomeron propagator now is a function of both rapidity and impact parameter. 
The pole contribution or the bare propagator reads 



It goes to delta function when the rapidity y goes to zero and its value at 6 = is 
growing with rapidity as ~ e^^/y. 

In figure 2 we show the calculation of the full propagator as a function of impact 
parameter in comparison with the bare one for different rapidities. The calculations 
show a negligible difference between sets of parameters which correspond to the same 
value of the triple coupling r^p and differ only by the value of the quartic coupling 
X ({1,2,3} and {4,5}). This supports the statement of [10] that presence of a finite 
number of higher than triple order terms in the RFT interaction Lagrangian may be 
ignored. 

For making comparison with zero-dimensional case we take the value of the prop- 
agator at 6 = multiplied by e = vra^ which in terms of the stochastic model is the 
average number of partons with position within that area and on the other hand is the 
amplitude for a single-Pomeron exchange. One may expect saturation in this quantity 
at sufficiently high rapidities even in presence of diffusion, since, naively speaking, at 
low densities parton splitting will prevail and at high densities fusion and annihilation 
will increase resulting in saturation eventually. However, as one can see from fig. 3 
the saturation is not reached within the range of rapidity achievable by the numerical 




(27) 
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Figure 2: The full Pomeron propagator as a function of impact parameter for different values of 
rapidity y in comparison with the perturbative one. 

calculation, unlike the case of zero transverse dimensions which is equivalent to zero 
slope. This is a consequence of the fact that the value of the slope which we used and 
which is close to the value dictated by the fits to the experimental data (e.g [14]) is in 
fact quite large. 

The applicability of zero dimensional approach depends on the interrelation of the 
parameters of the Reggeon Field Theory. Within the stochastic model the applicability 
condition can be formulated as follows. The number of splittings/fusions of partons 
must be already large for the rapidity ya which is required for the partons to diffuse 
to the distance of about the parton size. 

That is, within the stochastic approach the regime similar to zero- dimensional case 
is realized for 



This situation is illustrated in fig. 4 where we draw the same value as in fig. 3 for the 
parameter set number 1. However this time the curves correspond to different values 
of the slope a' ranging from 0.01 fm^ to 10~^ fm^ (from 0.26 to 2.6 ■ 10"^ GeV~^). 
One can easily see that the zero-dimensional regime is observed for a' = 10"^ fm^ up 
to ?/ ^ 10 which is in agreement with the estimate (28). However, since the triple 
Pomeron coupling is comparatively small (and, consequently, A and a are), it is rather 
the opposite limit that must be adopted for the applications and the diffusion must be 
fully taken into account. 

To conclude, in presence of diffusion the full Pomeron propagator at a fixed distance 
in the transverse space for the realistic slope value dictated by fits to the experimen- 
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Figure 4: Full propagator scaled by e = vra^ at 6 = for the constants A, u, mi and m2 of the set 
1 and different slope a' values. 
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tal data is a growing function of rapidity within the reach of numerical calculations. 
Within this rapidity range it just expresses a slower growth than the bare one (pole 
contribution). This is unlike the case with the absence of diffusion where the full 
propagator drops at high energies or reaches constant value comparatively fast. We 
also note that the behaviour of the propagator at large rapidities seem to depend on 
the slope value in a highly nontrivial way. This points out an important role of the 
Pomeron slope for the set on of saturation and probable change to a subcritical be- 
havior in the propagator. Another important distinction between the zero- dimensional 
and 2-dimensional theory is the role played by the quartic coupling. As one can see, 
the quartic coupling drastically changes the propagator asymptotic behaviour in case 
of zero-dimensions, while in presence of diffusion its role is negligible. 



5 Total and elastic cross sections. 

5.1 Normalization of the amplitudes and the couplings. 

To proceed with applying the technique to the calculations of the cross sections we first 
fix the normalization of the amplitude. This procedure is essential for the applications 
as it affects definitions of the phenomeno logical couplings of the Pomeron to particle. 
We use the elastic scattering amplitude M{Y, q) normalized in such a way that 



a'°\Y) = 2ImM(F,q = 0), 




where Y = ln(s/so), Sq = 1 GeV^ and q^ = —t. The amplitude in the coordinate 
representation f(Y, b) is defined as 

f{Y, b) = ^ J d^qe-'^^MiY, q) . (30) 

With this definition of the Fourier transform, the perturbative Pomeron propagator 
in the coordinate representation coincides with the Green function of the diffusion 
equation with splitting. The cross sections are 

a*°*(F) = 2 J d%lmf{Y,h) , = ^ d^lfiY,}))]"^. (31) 

We neglect the small real part of the Pomeron signature factor, so in what follows the 
Pomeron amplitudes are considered as purely imaginary, f{Y, b) ~ iT(Y, b), T = Im /. 

To establish the correspondence between the Pomeron-particle couplings and the 
distributions of projectile- and target-associated partons in the stochastic approach we 
use the contribution of single Pomeron exchange to the total cross section (fig. 5a). In 
the standard Regge parametrization the Regge-pole contribution has the form 

Tp(y, q') = (3^{q')(3^{q') e-"^«'^e^^ , (32) 

where are the couplings of the Pomeron to the colliding hadrons. They are often 
written in the form /3(g^) = f3og{q) with g{q = 0) = 1. The most common choice for the 
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(a) (b) (c) 

Figure 5: Reggeon graphs used for normahzation of the couphngs. 

profile g is the Gaussian: g{q) = exp(— In the impact parameter representation 
one has 



Tp{Y,b) = (3^(3^ 

j d%d% fi^^(bi) G{Y, b2 - h,)g^{h - b2) , 
where g^^^/x Fourier images of the couplings profiles normalized to unity 

fi^A/A(b) = j d'qe~'^''g^^x{q) , 

and 

62 N 



(33) 



(34) 



exp^ 



is the Pomeron propagator exp((A — «'q^)y) in the coordinate representation. 

To model a single Pomeron exchange in the stochastic model approach we set 
7722 = 1^ = and take only the term with s = 1 from the sum (15) (this corresponds 
to = 1 terms in (18)). For the imaginary part of the amplitude using (9) and (23) 
we evidently get 

Tp{Y,b) = eN^N^ 

J d%,d% p/hi) G{Y, b2 - bi)p^(b - ba) , (36) 

where Nj^^j stand for the average number of partons attributed to the projectile A and 

target A at the start of evolution and Pj^^/xi^) their distributions in the transverse 
plane which are normalized to unity. The function G{Y] b) describes diffusion of 
partons in the transverse plane and growth of their number with rapidity. It coincides 
exactly with the Pomeron pole propagator (35) if the diffusion coefficient D = a'. It 
follows from comparison of eqs.(33) and (36) that the initial parton distributions in the 
transverse plane coincide with the Pomeron-particle coupling profile in the coordinate 
representation and the number of partons is proportional to the value of the coupling 
at zero momentum transfer. 

(^o^%Aib) = VeN^/APA/xib) • (37) 
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The total cross section for the proton-proton case where A'^ = A^^ = is corre- 
spondingly 

= 2/3^6^3'. (38) 

where the constant /3 = /3q = N^/e. 

Now it is straightforward to write projectile- and target-associated parton distribu- 
tions in the quasieikonal approximation. The quasieikonal amplitude is written as an 
expansion in the number s of Pomeron exchanges with a special choice for the particle- 
s-Pomeron vertices. In the normalization of the amplitude which we have adopted it 
reads: 

T(r, h) = Y^ t^\y, h) = Y, — my. w ■ m 

s=l s=l 

The eikonal approximation corresponds to the special choice of C = 1. Writing down 
the s Pomeron exchange contribution for F — )■ explicitly and taking into account (9), 
one finds that it corresponds to the convolution of two s-parton inclusive distribution 
of the following form: 

e^/'ff/\h,, . . . , b,) = tlPA/Ai^^). (40) 

i=l 

This corresponds to the partons independently distributed in the transverse plane with 
the distribution Va i(bj) and parton number distribution with the factorial moments 

^ {,rr c'-^'^/U tL set of the factual .o,„e,.s its turn 

to the "quasipoissonian" distribution (poissonian distribution for the eikonal) with 

P;v = «:^-'^e-^^ n = l,...,oo (41) 
^0 = 1- E^=i = 1 - ^(1 - e"^^) , (42) 

where k = \/C and 7 = — 

Ve 

It is interesting to note that the quasieikonal approximation (39) can be reproduced 
within the stochastic approach for any rapidity. This is provided if the relation between 
the inclusive distributions, fsiv, bi, . . . , b^) = C^'^"^^^'^ Yl fiiVi bi), which valid at y = 
due to the independent distribution of partons in the transverse plane is also fulfilled 
for any rapidity value. This amounts to the factorial moments of the parton number 
distribution growing as 

/i,(l/) = e%,(y = 0) (43) 

and (for the Gaussian initial profile) independent Gaussian distribution of parton co- 
ordinates in the transverse plane with standard deviation growing along with rapidity 

as 

R\y) = R' + a'y. (44) 

To fulfill these conditions in the course of computation it is necessary not only to 
forbid fusion and annihilation of partons but to impose an additional condition on 
the level of linking of the inclusive distributions fs{Zs). Namely, all partons which 
coordinates enter the set Zs must belong to different connected components, that 
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Figure 6: Schematic view of the contribution entering the quasieikonal amphtude in the stochastic 
approach. Only graphs of type c) contribute. 

is, must be produced from different initial partons. In our numerical scheme this 
amounts to taking only such contributions to the sum in the (12) which do not contain 
coordinates of the partons from the same connected component. 

The corresponding numerical procedure is realized by dividing the matrix G = \\gik\\ 
entering (13) into blocks which correspond to interaction of partons belonging to the 
same connected components. The blocks are numbered by connected components of 
the graph which come from = and y = Y . Terms entering the sum in (13) are 
those whose multipliers come from different blocks with different components number, 
each number of component can not be used twice in same product. In other words, 
all the multipliers must come from different block lines and different block columns. 
Effectively this corresponds to substituting the matrix G by G which dimensions are 
the numbers of connected components and elements being sum of the elements of the 
matrix G which enter the corresponding block. The matrix G is processed according 
to the usual procedure described in the section 3. 

The contributions which enter the amplitude (13) computed on the partonic samples 
can be schematically illustrated by the fig. 6 where solid lines correspond to parton 
propagation in time. Obviously, in the quasieikonal approximation we sum up only 
a certain part of the contributions to the amplitude. However the result is obviously 
independent from the position of the partonic samples linkage point in rapidity, the 
same as for the full calculation with account of fusion and annihilation. 

As the next step we establish a link of the triple Reggeon coupling, r^p = \y/e in 
our notations, with the relevant quantities used in fits to the experimental data by 
different groups (see [11], [12], [13]). To do this we write explicitly the expression for 
the diffractive cross section in the lowest order in parton splitting, which is expressed 
via the single triple interaction contribution to the amplitude (see fig. 5c). Within the 
stochastic model we start with the contribution to the amplitude in the Schwimmer 
approximation with a' ^ which is first order in r^p. The Schwimmer approximation 
implies eikonal vertices for the target and m2 = 1^ = 0. In this case it is straightforward 
to write down a solution to the set of equations (43) of [3] for the inclusive distribu- 
tion. The convolution of f2{y, bi, 62) with the particle - two Pomeron vertex gives the 
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stochastic model expression for the amphtude in the impact parameter representation: 
Tsp{Y,b)= r dyfsp{y,Y;b) = 

^ 2/min 

-Xe^N^Nl r dy ! d%Ad%2dXpA{hi)G{yM-hi) 

X G{Y -yM- bo)G'(F - - 6o)l>i(b - 62K(b - fo'a) . (45) 

Its discontinuity which corresponds to the intermediate state with rapidity gap y 
(diffractive dissociation) is of the positive sign and double magnitude. The corre- 
sponding single diffractive cross section is 



d^ ...da' 
dyM 



dcT^^ r 
M'^ = -2 j d%T,p{yM,Y;b) 



= r3./3^(0)e^(^-^) J ^ (/3^(t))^^(^+"'*)^ . (46) 

where yM = Y — y = ln(M^/so). 

We make direct comparison of (38) and (46) with the corresponding definitions of 
[11] (rx), [12] (re) and [13] (r^) to find out that 

^ [T re 

rap = \/T^rK = \hr^L 



27r 2^2 ' 
5.2 The cross sections. 

In order to estimate the role of loop corrections it is instructive to compare our cal- 
culations with one of the realistic fit to the experimental data. Our aim is not to fit 
the data better - we want to estimate the relative contribution of the loop diagrams 
at realistic parameter values. As a starting point for the numerical calculation of the 
cross section we take a quasieikonal fit to data according to [14] with 

T(y, b) = ^ ^ {Tp{Y, h)r , (47) 



n=l 



where 



Tp(r, b) = ^^1^ exp[-l6V(i?l + o.'Y)]. (4J 



The fit has a Gaussian parametrization of the Pomeron-particle vertex. The parame- 
ters of the fit are as follows: 



A 


= 0.12 , 




9l 


= 2.14 GeV-2 ^ 


0.083 fm^ , 


a'p 


= 0.22 GeV"^ R 


i 0.0085 fm^ 


Rp 


= 3.30 GeV^^ ^ 


0.128 fm^ , 


C 


= 1.5 . 





(49) 



18 



We take initial parton distribution in the transverse plane according to this fit in the 
Gaussian form accounting for a different amplitude normalization in [14] ^: 

and the factorial moments of the parton number distribution as 

^ ^ a 

We fix the Pomeron intercept A = A — mi = 0.12 and the triple Pomeron coupling 
= Xy/Tw? = 0.087 GeV~^ which corresponds to the value = 0.05 GeV~^ from 
[11] in our normalization. 

Our aim is to illustrate dependencies of the cross sections on the parton size (linkage 
distance) a and quartic coupling x- Following the relations outlined in table 1 and 
throughout the section we define several sets of the partonic scheme parameters which 
correspond to the parameter set (49) and differ with transverse size of the parton a 
and the values of the quartic coupling x- We fix the parton size values at the level 
of 5% and 10% of the proton radius. The quartic couplings depending on relation 
between A and m2 for the sets are Xi = X4 = 0.0005569 fm^, X2 = 0.0002785 fm^, 
X3 = 0.0011134 fm^. The sets of parameters are presented in table 3. 

Table 3: Parameters used for the full calculation within the stochastic scheme. For all the sets 
R = 0.36 fm and D = 0.008538 fm^ are taken. 



Set 


a, fm 


A 


mi 


m2 


V 




1 


0.018 


0.54722 


0.42722 





1.09488 


32.02 


2 


0.018 


0.54722 


0.42722 


0.54722 





32.02 


3 


0.036 


0.27361 


0.15361 





0.54722 


16.01 


4 


0.036 


0.27361 


0.15361 


0.27361 





16.01 



The results of the full calculation are presented in figure 7. As one can see, Pomeron 
interaction leads to a certain shadowing in the amplitudes and hence slows down the 
cross section growth compared to the quasieikonal approximation. The value of the 
cross section is mainly defined by the regularization scale (parton size a) and not by the 
value of the quartic coupling, the same as for the Pomeron propagator. We have also 
checked numerically the independence of the result on the linkage point and found it to 
hold within 2% accuracy. We have also verified that the quasieikonal approximation is 
reproduced numerically in our computation once we impose restrictions on the parton 
sets evolution and their linkage procedure as described in section 5.1. 

It is possible by changing parameters to bring the results of the full calculation 
close to the quasieikonal fit to data. We did so by leaving all the coupling the same 
and increasing the intercept value simply reducing mi. As one can see in the fig. 7 
changing the value of A = A — mi to 0.165 brings the full calculation results close to the 
original quasieikonal fit (A = 0.12). The curves for the derived sets with A = 0.165 are 

^In paper [14] one uses normalization (t*°*(s) = 87rlmr(s, q = 0). 
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additionally marked with an asterix (1*, 2*, 3*, 4*) with respect to the original sets 
with A = 0.12. The change in the value of the quartic coupling x has a minor influence 
on the cross section, much less than a change in the regularization scale (parton size) 
a. 

The same as for the intercept value 1.12, the quartic coupling value has a minor 
influence on the values of the cross sections. For the energies reachable in the nearest 
future changing quartic coupling by a factor of 2 amounts to less than 5% change in 
the cross sections. 

6 Conclusion. 

We developed a numerical realization of the RFT which allows to obtain the exact value 
of scattering amplitude taking all loops into account using the equivalence between 
the Reggeon Field Theory with Pomeron scattering term and parton stochastic model 
(reaction-diffusion approach). Unlike other calculations with account of loops for no- 
zero Pomeron slope [13, 15, 16] our aproach does not involve inflnite number of flne- 
tuned couplings for transitions of m into n Pomerons where results depend strongly 
on the chosen ra, m dependencies of the vertices. Our calculations includes triple and 
quartic couplings only and, to our knowledge, is the flrst numerical all-loop calculation 
for the case of non zero slope. This calculation possesses the expected properties of 
the approach: boost invariance and projectile-target symmetry. Moreover, it gives 
the exact amplitude as an expansion in the number of the Pomeron exchanged at a 
given rapidity value y which provides an straightforward opportunity of calculating 
the quantities which involve this expansion, such as multiplicity distribution in the 
narrow rapidity interval. 

The stochastic model has several parameters, most of them have direct correspon- 
dence within RFT and are flxed through its couplings. One of the two parameters 
which are not flxed at once is the parton interaction distance a. It is equivalent to 
the ultraviolet cutoff or the Pomeron size for the calculations within the RFT thus 
playing a role of the RFT renormalization scale. Another parameter is the ratio of 
parton fusion to parton annihilation probabilities. This ratio deflnes the relative value 
of triple and quartic Pomeron couplings and up to our knowledge there are no ex- 
perimental constraints on its value yet. However, as we show numerically the role of 
the quartic coupling in the RFT with the slope a' distinct from zero is negligible for 
the asymptotic behaviour of the amplitudes and the cross sections within the reach of 
numerical computations. This is unlike the case of the zero slope. 

As an illustration of use of our approach we compute the total and elastic proton- 
proton cross sections in the approximation of quasieikonal proton-Pomeron vertices. 
Comparison of the full calculation with the quasieikonal flt shows that for the relevant 
energy range the full account of the Pomeron interactions can be effectively imitated in 
the quasieikonal approximation by reducing the intercept from A = 0.165 to A = 0.12. 

We hope that our study will be helpful in understanding the role of relative value 
of the Pomeron loop contribution in high energy hadronic scattering. 
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Figure 7: Total and elastic cross sections; full calculation compared to the quasieikonal fit [14]. 
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